7T7T Scattering in Twisted Mass Chiral Perturbation Theory 
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In this report we describe both 1 = 2 and 1 = nn scattering for twisted mass lattice QCD 
utilizing twisted mass chiral perturbation theory at next-to-leading order. Focusing on the lattice 
spacing (b) corrections, we demonstrate that in the exotic / = 2, 13 = ±2 channels (jr n ), the 
leading scaling violations of -kit scattering at maximal twist begin at 0(m%b 2 ). This is not the 
case in any other isospin channel, for which the scaling violations at maximal twist begin at 0(b 2 ). 
Furthermore, we demonstrate the existence of a mixing between the I = 2, 13 = and 1 = 
scattering channels due to the breaking of isospin symmetry by the twisted mass term. The mixing 
term, although formally next-to-leading order, is relatively large, thus necessitating the use of a 
coupled channel analysis. We argue that this mixing likely renders the computation of the 1 = 
channel impractical with twisted mass lattice QCD. 

PACS numbers: 12.38.Gc 
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INTRODUCTION 



The last few years have seen a growth in the computa- 
tion of hadron interactions with lattice QCD fl], with ( 
namical calculations of two- meson systems 0, 1, El H, 
two-baryon systems @, [H and systems of up to 12 pi- 
ons [H, |T(| and kaons Further, lattice field the- 

ory methods are now being applied to the low-energy 
effective field theoryof multinucleon interactions [12], [l3|, 
HI EEES d El, tH, for which there exists a nice re- 
view [2(|. The dynamical lattice QCD calculations of 
hadron interactions to date, have either been performed 
with Wilson fermions or a mixed lattice action [U [22j of 
domain-wall valence fermions 0, [24], [25[ and the Asq- 
tad improved [26l . \2l\ rooted staggered MILC configura- 
tions [SHH. Twisted mass lattice QCD [HSU has re- 
cently emerged as a viable fermion discretization method 
for computing gauge configurations with two flavors of 
light quarks, (up and down) [H, [H, HI] and hopeful 
prospects of 2 + 1 + 1 (up, down, strange and charm) 
flavors of dynamical sea fermions (35| in the chiral regime. 
It is therefore only a matter of time before hadron inter- 
actions will be computed with the twisted mass fermion 
discretization method. 

The first multihadron system to be explored with 
twisted mass lattice QCD will most likely be that of two 
pions. The two-pion system is numerically the simplest 
as well as theoretically the best understood. In fact, the 
scattering of two pions at low energies was uniquely pre- 
dicted at leading order (LO) in chiral perturbation theory 



(%PT) by Weinberg in 1966 36]. The subleading orders 
in the chiral expansion give rise to perturbative correc- 
tions to the LO predictions and have been worked out to 
one-loop, or next-to-leading order (NLO) by Gasser and 
Leutwyler [3 71 ] and also to two- loops, or next-to-next- 
to-leading order (NNLO) [H, [H Sfl]. Each new order 
introduces operators with coefficients not constrained by 
chiral symmetry, known as low-energy constants (LECs). 
To have predictive power, these LECs must be deter- 
mined either by comparison with experiment or lattice 
QCD calculational results. A comparison with lattice 
QCD can introduce additional complications as the cal- 
culations are performed at finite lattice spacing in a fi- 
nite volume. Modifications to the infrared and ultravio- 
let behavior of the theory can be incorporated into chiral 
perturbation theory. For sufficiently large but finite lat- 
tice volumes, the operator structure and power counting 
of the effective theory remain valid with exponentially 
small corrections to matrix elements [4l|. Lattice dis- 
cretization effects can also be incorporated into the chi- 
ral Lagrangian through a two-step process first detailed 
in Ref . 0] ■ One first constructs the effective continuum 
Symanzik Lagrangian [43|, [44[ for a given lattice action. 
One then builds the low energy chiral Lagrangian from 
the Symanzik theory, giving rise to new unphysical oper- 
ators with their own LECs. These new operators capture 
the discretization effects for a given lattice action. 
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In this report, we briefly review the construction of the 
twisted mass chiral Lagrangian in Sec. [Til We then de- 
termine the lattice spacing (b) corrections to low-energy 
7T7T scattering specific to the twisted mass lattice action. 
We work through 0(b 2 ) - 0(bml ). 



2 



II. TWISTED MASS LATTICE QCD AND THE 
CONTINUUM EFFECTIVE ACTION 



The twisted mass chiral Lagrangian was determined 
previously in Refs. [H, IH, E3, S3, H§| , and for baryons in 
Ref. [5(| • In this report we focus on twisted mass lattice 
QCD with degenerate light flavors given by the lattice 
action 

s = E fa) \l E >( v ^ + v -) - r o E v ^ 



+ to + «75T3Mo 



tp(x) 



(1) 



where ij) and tp are the dimensionless lattice fermion 
fields, V„(V*) are the covariant forward (backward) lat- 
tice derivatives in the v direction, too is the dimensionless 
bare quark mass and /io is the dimensionless bare twisted 
quark mass. The fermion fields are flavor doublets, t 3 is 
the third Pauli-spin matrix and the bare mass term is 
implicitly accompanied by a flavor identity matrix. Our 
twisted mass xPT analysis also holds for dynamical lat- 
tice calculations with 2+1 + 1 flavors, the only differ- 
ence being the numerical values of the LECs determined 
when fitting the extrapolation formula to the calculation 
results. 



The continuum chiral Lagrangian, supplemented by 
discretization effects is determined with the two step pro- 
cedure of Ref. 1421. This was done for the twisted mass 



lattice action in Ref. |47j , to NLO in which a power count- 
ing vriq ~ bAQ CD was used and which we shall adopt. 
The resulting effective Lagrangian is 

C eff = Cgi ue + q{Ip+ to + ij 5 T 3 /j,)q 

+ c S wb qia-^F^q + 0(b 2 ,bm q ,m 2 q ) , (2) 

where C g i ue is the Yang-Mills Lagrangian. The quark 
fields are an isodoublet, q T — (q u ,qd) and the quark 
masses are given by 



m = Z m (m -m c )/b, 



(3) 



The symmetry properties of the twisted mass lattice ac- 
tion protect the twisted mass from additive mass renor- 
malization. With Eq. ([5]), one can construct the two 
flavor chiral Lagrangian. This is the Gasser-Leutwyler 
Lagrangian (37j supplemented by chiral and flavor sym- 
metry breaking terms proportional to the lattice spac- 
ing. The Lagrangian through NLO relevant to our 
work takes the form [13, |49[ (we use the normalization 
/~130MeV), 



4™ = ^MfyE^Et) - ^-tr( X 't£ + EV) - ^(cVEc^Et) 2 - ^tr^EctE^tr^Ec^Et) 



r- 



h. 



h 



h + k 
16 



tr(x /f E + SV) 



k. 



tr(9 M E9 M Et)tr(x' t E + EV) 
VKtr(^ I Ea M E t )tr(i t E + E f i) - Wtr(x' f E + EV^r^E + E f i) - W [tr(i f E + EU) 

I 



(4) 



where the LECs, I1-I4 are the SU{2) Gasser-Leutwyler 
coefficients and the coefficients W, W and W' are un- 
physical LECs arising from the explicit chiral symmetry 
breaking of the twisted mass lattice action. The spurion 
fields are defined as 



X 1 = 2B Q (m + ir 3 n) + 2W b = m 
A = 2W a b= b. 



(5) 



As discussed in Ref. [49(, the vacuum of the theory as 
written is not aligned with the flavor identity but is given 
at LO by 

E = (0|E|0) = — = exp(iw r 3 ) , (6) 

with 



Therefore, to determine the Feynman rules which leave 
the interactions of the theory the most transparent, one 
expands the Lagrangian around the physical vacuum. 
Extending this analysis to NLO, one finds the vacuum 
angle shifts to lo = luq + e where one can determine e 
either by finding the minimum of the potential, as was 
done in Ref. [491 ] or by requiring the single pion vertices 
to vanish, 



32 - 

e(w ) = --^fcsinwo 
f 



W + 2W' cos u — 



(8) 



One can expand about the physical vacuum by making 
the replacement 



M' = 



\J \rh + b) 2 + /i 2 . 



(7) 



E = £ m E p/l £ m , with £ m = exp(iwr 3 /2), (9) 
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and 



One then finds the Lagrangian is given by 



T, ph = exp [ — 




(10) 



73/ 



C = £ c 



WbcosLjti(d^ ph d^l h )ti(J: ph + Ejj - fccosw (Vm' + W'Scobw) [tr(E p/l + E+J 

wSBinwtr^Ep/ifl^EjfcJtrCiTaCSpfc - E^)) - W6 2 sin 2 w [tr(iT 3 (S p/l - Ej fc )) 

M'/ 2 



tr(ir 3 (S p /, - E^)) 



e( W )- 



6 sin w (WM 1 + 2W r bcosuj^ tr(E p/l + T, j ph ) 



(11) 



r 



where C CO nt. is the continuum (2) chiral Lagrangian 
to NLO. Of particular interest to us are the new two, 
three and four pion interactions which result from the 
discretization errors in the twisted mass Lagrangian. We 
find, in agreement with Ref. [4^ | 



C=C B 



A£ 24> + AC 3<b + AC 



30 



-40 : 



(12) 



where 



AC 2l p = cosw^^tr^d^) + ]-AM'{uj)tv{4> 2 ) 
J 2 



tr 



T30 

V2 



n 2 



(13) 



16Wb 
1 

e(uj)M' 



A£ 3t p = - sinw tr(T 3 0)tr(9 M< /)g At (/)) 



2/ 



tr(r 3 ^)tr(^) . 



(14) 



AM>) 2 2 
A£ 40 = — — [tr(0 )J 



3/ 2 



3/ 2 



tr(^)] 2 tr(0 2 ) 
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16Wb , n ,„ ,. . , 9 . 
-coscj tr(ej M 0d M 0)tr(0 ). 

and the mass corrections are given by 



(15) 



Al/'u) = cnsw^- ^M' + coswlf'6 
64iy'S 2 



AM' Q {u) = -sin 2 w- 



P 



(16) 



From this Lagrangian, one can determine the pion 
masses, decay constants and wave-function corrections. 



One finds the masses arc (using the modified dimensional 
regularization of Ref. (37j ) 



m 2 ± = M' 



1 



M' , /M'\ „. ,4M' 



(4tt/)- \ /x 

32W&M' 



+ AM'(w) - cosw 
m 2 = m 2 ± + AM' Q {w) , 
the decay constants are 1 



/ 2 



(17) 
(18) 



2M' , /M'\ ir/ ,2M' 



(47T/)* 



/ 2 



16W^6 



/ 2 



and the wave-function correction is 



(19) 



AM 1 , / M' 
In 



3(4*/)* 



Km) 



4M' 



32VK6 



/ 2 



(20) 



These expressions will be needed to express the scatter- 
ing in terms of the lattice-physical parameters (by lattice- 
physical, we mean the renormalized mass and decay con- 
stant as measured from the correlation functions, and 
not extrapolated to the continuum or infinite volume 
limit). As we discuss in the next section, these inter- 
actions lead to three types of new contributions to 7T7t 
scattering states: there are discretization corrections to 
the scattering parameters, the scattering lengths, effec- 
tive ranges, etc., which appear in a mild manner as those 
from the Wilson chiral Lagrangian [H3|. There are cor- 
rections which can potentially significantly modify the 



1 There is an exact Ward identity one can exploit to compute the 
charged pion decay constant and avoid issues of the axial current 
renormalization discussed for example in Ref. fsitl . 
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(a) 



(b) 



FIG. 1: New unphysical graphs from twisted mass interac- 
tions in the £(it)-channel (a) and s-channel (6). Fig. (6) can 
only contribute to /3 = scattering. 



chiral behavior, arising from the three-pion interactions, 
and there are new corrections which mix different scat- 
tering channels, for example, the / = 2, I3 = and the 
1 = scattering states. 



III. 7T7T SCATTERING IN TWISTED MASS xPT 

In this section we calculate corrections to the two- 
pion scattering channels. We begin with the maximally 
stretched 1 = 2 states, which have the simplest correc- 
tions. 



A. I = 2,h = ±2 Channels 

There are two types of discretization corrections which 
modify the / = 2, ^3 = ±2 scattering, those which 
are similar to the corrections for the Wilson lattice 
action [HI, [HHJ and those which arise from the three- 
pion interactions, Eq. (fT4f and give rise to new Feyn- 
man diagrams. We will express the scattering parame- 
ters in terms of the lattice-physical pion mass and de- 
cay constant. As was shown in detail, this has dra- 
matic consequences on the formula for the scattering pa- 
rameters in both partially quenched and mixed action 
xPT [H(| [57], HH], such that the extrapolation formu- 
lae were free of unphysical counterterms through NLO. 
There is a second benefit to expressing the scattering 
parameters in lattice-physical parameters. This allows 
one to perform a chiral extrapolation in terms of the 
ratio m^/f^, and thus avoid the need for scale setting. 
This was crucial in allowing the NPLQCD Collaboration 
to make a precision prediction of the 1 = 2 scattering 
length @,[|. 

The simple corrections to the scattering amplitude are 
determined from AC^, Eq. (fT5)) . The three-pion inter- 
actions from Eq. (| 14[) lead to new topological graphs in 
the scattering amplitude, which we depict in Fig. [TJ The 
I = 2.I3 = ±2 scattering channels receive corrections 
from Fig. [IJa) and its u-channel counterpart. The inter- 
nal propagating pion is a 7r°, which, for present lattice 
actions, is known to be lighter than the charged pions in 
dynamical twisted mass lattice calculations with degen- 
erate light quark masses [HI, [34[ . 

Putting all the corrections together, one finds the scat- 
tering amplitude, which we express in Minkowski-space, 



is given by 



n-2,±2 _ t-2,±2 A ^-2,±2 
1 — 2 cont. ' 



(21) 



where the discretization corrections are 



AT 2 ' ±2 (u) = 4 



AM'(w) 



cos ui — — — s 



f 4 



2t 2 {uj)m% 

J TV 



t 



mlo - u 



(22) 



The first two terms arise from Eq. (fTSj) as well as from 
the conversion to the lattice-physical parameters. The 
second two terms arise from Fig. (Ua). These terms are 
formally NNLO. However, depending upon the precision 
with which the twist angle is tuned, these terms may 
become large and require promotion to lower order. Ex- 
panding the NLO contribution to the twist angle, Eq. ([3]), 
one finds 



64W'b 2 



• 2 2 

sin to cos to 



+ 0(ml). (23) 



We can then determine the corrections to the I = 2, 13 = 
±2 scattering lengths, for which we find 



A? 



,1=2, ±2 



M = 



AM'(w) 
^7T 



— cos(w) 



SWbml 



+ 



(32FT') 2 sin^cos 2 ^ b 4 



(24) 



The first observation we make is that at maximal twist, 
lu = 7r/2, these leading discretization errors exactly can- 
cel through NLO (this is true of the corrections to the 
scattering amplitude and not just the scattering length) 2 



Am x a^ ±2 (V2) = 0. 



(25) 



This is independent of the use of lattice-physical parame- 
ters, and holds also for the scattering length expressed in 
bare parameters, or any combination of bare and physi- 
cal. At zero twist, u> = 0, our expressions reduce to those 
of Ref. [H3]. Converting f n — > /, our answer agrees with 
that in Ref. [HH| . The scattering length at maximal twist 
is simply given by the continuum formula 



1=2 



-2tt 



f m ^ 



3 In 



(26) 



2 We have assumed that a suitable definition of the maximal twist 
angle has been used in the numerical lattice computations such 
that one is not restricted to the regime m q >> b 2 A^ CD , but 

El, El nim. 



rather one is allowed m q > 6Aq CD 
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where the combination of Gasser-Leutwyler coefficients 
is Sill 



In the isospin limit, the |2,0) and 1 0, 0) states (in the 
\I, I3) basis) are given by 



l'= 2 = 4(^(411 + Al r 2 +l r 3 -Q. 



(27) 



Furthermore, the discretization errors only enter at tree 
level at this order (when the expression is expressed in 
lattice-physical parameters), and thus at arbitrary twist, 
the exponentially suppressed finite volume corrections to 
Liischer's method are also given by those determined in 
continuum finite volume x^T (59j . 

Returning to the new graphs arising from the three- 
pion interactions, we can estimate the size of the correc- 
tions to the scattering amplitude using the known mass 
splitting between the charged and neutral pions [H, HH • 
Estimating the splitting with the leading correction, 
Eq. (fT5|) , and solving for W' from Eq. (|16[) . we can es- 
timate the corrections to the / = 2,1$ = ±2 scattering 
length near maximal twist. As a ratio to the LO predic- 
tion for the scattering length, one finds 



Am w ai= 2 ' ±2 M cot 2 u (AM^)/m 2 



,7=2, ±2( 



l + AM^(w)/m| 



(28) 



At the lightest mass point calculated in Refs. [32|, l34j |. 
which corresponds to m w ~ 300 MeV, the pion mass 
splitting is 



AM^u ~ tt/2) 



-0.33. 



(29) 



and therefore one must have cotu > 0.3 for this term 
to make more than a 1% correction. Therefore, for cur- 
rent twisted mass lattice calculations, corrections to the 
I = 2, 13 = ±2 scattering length (and other parameters) 
should be negligible provided higher order corrections are 
as small as expected. 



B. J3 = scattering channels 

There are several features which make scattering in 
the J3 = channels more complicated than in the 
I = 2 ; Is = ±2 channels, most of which stem from the 
fact that the twisted mass lattice action explicitly breaks 
the full SU(2) symmetry down to t/(l), the conserved 
I3 symmetry. The first technical complication is not 
specific to twisted mass calculations, but is simply the 
need to compute quark disconnected diagrams. The sec- 
ond complication stems from the mass splitting of the 
charged and neutral pions. Generally, one determines the 
scattering phase shift for two particles with the Luscher 
method [60, HH [si, [iiH , by determining the interaction 
energy 



|2,0) = ^=(k+^) + |7r-7r+)-2|7r% )) , 

|0, 0) = -L (K+TT-) + I7T-7T+} + IttV)) . (31) 

However, given the relatively large mass splitting in 
current twisted mass lattice calculations, Eq. (|2T)j). the 
propagating eigenstates will be arbitrarily shifted from 
the physical states, perhaps shifting nearly to the 
{|7r ± 7r T ), |7r°7r )} basis. This would have to be disen- 
tangled numerically. Even ignoring this issue, which we 
deem the most serious, and working with the continuum 
{|2, 0), |0, 0)} basis, there is a mixing of these states, 
which first appears at NLO as the second operator in 
Eq. (fT5|) . Working with the states 



\I,0) = 



2,0) 
0,0) 



one finds 



AT 20;0,0 



9/ 2 



4 -4s 

7 

~V2 



(32) 



(33) 



Given the correction AA/q(w), one sees this mixing is in 
fact maximal at maximal twist. This is nominally a NLO 
effect, thus possibly leaving the states mostly unmixed. 
However, a comparison of this term with the LO ampli- 
tude of the / = 2, J3 = ±2 scattering, one finds close to 
maximal twist 



ATi£ ;7(327r) 



1 



-27r(m w /4Tr/ T )2 TO 2 // 2 



1.11 -1.37 
-1.37 1.38 



(34) 



AE„ 



2V5 



2m„ 



(30) 



For m-nj fn = 2, all terms in this scattering matrix are 
approximately 1/3 the size of the LO amplitude. Since 
we now know that the physical NLO corrections to the 
1 = 2 scattering length for example, only provide a few 
percent deviation from the LO term 0, 0, [56|, , we 
conclude that this NLO operator in fact provides a rela- 
tively large contribution to the scattering amplitude, and 
furthermore provides a large mixing term, and thus can- 
not be neglected. This, combined with the problem we 
mentioned previously, means a coupled channel version of 
Liischer's method of determining the scattering param- 
eters would be needed to explore the 1% = scattering 
channels with twisted mass lattice QCD. 

The AC^cf, Lagrangian is not the only source of mix- 
ing. The three-pion interactions, depicted in Fig. [1] in 
the s, t and u channels, will also lead to a mixing of the 
|2,0) and |0,0) states, as one can check with an explicit 
calculation. One may be concerned that the new s chan- 
nel graph will invalidate Liischer's method. This is not 
the case however, as the internal pion propagator is al- 
ways off-shell, and thus these diagrams do not contribute 
to the power-law volume dependence of the two-particle 
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X 



(a) 



X + "W^ 7 + +... 
"^N. ^ 

(b) 

FIG. 2: Modified four-point function, (a) consisting of all off- 
shell graphs. These vertices can then be iterated and summed 
(b), to determine the irn interactions. This summation gives 
rise to the Liischer relation, valid below the inelastic thresh- 
old. 



energy levels. An alternative way to understand this is 
diagrammatically. One can define a modified (momen- 
tum dependent) four-point function, which is order by 
order all the diagrams which do not go on-shell below 
the inelastic threshold. We depict this modified vertex 
in Fig.[2ja). These 2PI diagrams can then be resummcd 
to all orders to produce the scattering matrix, Fig. &[b). 
It is this resummation that produces the Liischer rela- 
tion, relating the finite volume scattering to the infinite 
volume scattering parameters [10, HH H2, [Hj]. In this 
way, one can see that the new interactions will not lead 
to a modification of the structure of the Liischer relation. 

Our last note of caution regards the construction of 
the interpolating fields. In the physical basis, the states 
which become those of definite isospin in the continuum 
limit are given by Eq. (|3ip . However, the interpolat- 
ing fields are generally constructed with quark fields in 
the twisted basis, with a known definite twist from the 
physical basis fields. While this is also true of the tt + tt + 
scattering channel, the phase is trivial since there is only 
one term contributing to the |2,2) state. Thus if one 
were to undertake a calculation of these coupled scatter- 



ing channels, care should be taken in constructing the 
correct interpolating fields. 



IV. CONCLUSIONS 

In this report, we have detailed tttt interactions in 
twisted mass xPT. We have shown that through NLO, 
at maximal twist the corrections to the 1 = 2, 13 = ±2 
scattering parameters from discretization errors are iden- 
tically zero. However, near maximal twist there are cor- 
rections which can modify the expected chiral behavior 
which we demonstrated by an explicit calculation of the 
correction to the scattering length. We found however, 
that for the dynamical twisted mass lattice configurations 
which exist today, the expected corrections are negligible. 

The I3 = scattering channels proved to have more 
significant discretization corrections, most notably a mix- 
ing term between the |2, 0) and |0, 0) states which is rel- 
atively large. In fact, these mixing terms combined with 
the need for computing quark disconnected diagrams and 
the expected nonperturbative shift of the twisted mass 
eigenstates, as discussed in Sec. IIIIB1 may make a cal- 
culation of these -Z3 = scattering channels prohibitively 
complicated. 
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